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Abstract 

A  military  arms  race  is  characterized  by  an  iterative  development  of  measures  and  countermea¬ 
sures.  An  attacker  attempts  to  introduce  new  weapons  in  order  to  gain  some  advantage,  whereas 
a  defender  attempts  to  develop  countermeasures  that  can  mitigate  or  even  eliminate  the  effects  of 
the  weapons.  This  paper  addresses  the  defender’s  decision  problem:  given  limited  resources,  which 
countermeasures  should  be  developed  and  how  much  should  be  invested  in  their  development  so 
as  to  minimize  the  damage  caused  by  the  attacker’s  weapons  over  a  certain  time  horizon.  We  for¬ 
mulate  several  optimization  models,  corresponding  to  different  operational  settings,  as  constrained 
shortest  path  problems  and  variants  thereof.  We  then  demonstrate  the  potential  applicability  and 
robustness  of  this  approach  with  respect  to  various  scenarios. 
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1  Introduction 


The  term  arms  race  is  typically  used  to  describe  military  buildup  efforts  by  countries  that  are  in 
conflict  with  one  another  (e.g.,  US  vs.  the  USSR  during  the  cold  war  era,  India  vs.  Pakistan, 
Greece  vs.  Turkey,  etc).  Studies  of  such  phenomena  typically  address  strategic  issues  and  are 
commonly  found  in  the  political,  economics  and  strategic  planning  literature  -  see,  for  example, 
the  recent  study  of  military  technology  races  in  [12].  In  contrast,  the  present  paper  addresses 
operational  aspects  that  arise  in  arms  races.  In  particular,  it  focuses  on  how  much  resources  should 
be  invested  by  a  defender  in  such  a  race  and  how  to  time  these  investments.  We  consider  an  arms 
race  between  two  asymmetric  parties:  Red  (R)  and  Blue  (B).  R  is  the  attacker,  who  is  trying  to 
develop  an  assortment  of  new  weapons  to  attack  the  defender  B.  Being  aware  of  R7  s  capabilities, 
intentions  and  activities,  B  is  trying  to  develop  countermeasures  (CMs)  that  will  mitigate,  or  even 
neutralize,  the  effects  of  R7 s  weapons.  The  CMs  may  be  technological,  tactical,  or  both. 

If  R  completes  the  development  of  a  certain  weapon  and  makes  it  operational  before  B  is  ready 
with  appropriate  CMs,  then  R  inflicts  a  certain  damage  on  B  (typically  measured  in  casualties  and 
economic  damages)  per  each  time-unit  until  an  appropriate  CM  becomes  operational.  If  B  wins 
the  race  and  a  CM  is  operational  before  R  deploys  a  weapon,  then  the  damage  to  B  is  smaller 
when  that  weapon  becomes  available.  If  B7s  CMs  are  perfectly  effective  against  that  weapon  (see 
[5]),  the  damage  to  B  can  be  as  low  as  zero.  Given  a  set  of  existing  and  potential  weapons  to  be 
deployed  by  R,  the  problem  that  B  faces  is  how  to  utilize  its  limited  resources  to  develop  the  most 
effective  mix  of  CMs  -  a  mix  that  minimizes  total  damage. 

An  example  of  the  settings  addressed  in  this  paper  is  the  counterinsurgency  warfare  faced 
by  coalition  forces  in  Iraq  and  Afghanistan  (2003-2009)  where  the  insurgents  develop  and  deploy 
new  types  of  improvised  explosive  devices  (IED),  with  ever  increasing  lethal  capability,  while  the 
coalition  forces  continue  to  develop  technologies,  tactics,  techniques  and  procedures  to  respond  to 
that  threat  (see,  e.g.,  [14],  [16],  [17]). 

Arms  race  problems  are  related  to  a  broader  class  of  problems  addressing  investment  rates 
in  R&D  projects  that  are  carried  out  in  competitive  market  environments  (see,  e.g.,  [9],  [11]  and 
[15]).  Most  of  the  articles  that  have  appeared  in  this  literature  have  assumed  the  “winner-takes-all” 
hypothesis  whereby  the  first  party  that  achieves  an  advantage  maintains  it  indefinitely  and  all  other 
parties  lose. 

In  a  recent  paper,  Golany  et  al.  [8]  analyze  a  stochastic  version  of  the  arms  race  problem  of 
the  kind  described  above.  In  contrast  with  the  common  “winner-takes-all”  assumption,  the  models 
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presented  in  [8]  address  situations  in  which  any  advantage  gained  by  one  of  the  parties  partici¬ 
pating  in  the  race  is  temporary  in  nature  and  is  lost  once  another  party  overtakes  the  lead.  Two 
types  of  models  are  presented  in  [8]:  optimization  models  that  derive  optimal  resource  allocation 
schemes,  and  game-theoretic  models  that  derive  N ash-equilibrium  solutions.  Specifically,  resources 
invested  by  B  in  developing  CMs  determine,  probabilistically,  the  time  when  these  CMs  are  ready 
and  operational,  and  thus  also  determine  the  expected  damage  caused  to  B  by  R.  Assuming  a 
predetermined  development  policy  of  CMs  -  in  parallel  or  sequentially  -  corresponding  convex 
programming  problems  are  formulated  and  solution  methods  are  discussed. 

The  stochastic  models  in  [8]  enable  B  to  determine  optimal  investment  schemes  while  capturing 
uncertain  durations  of  R&D  activities  and  limited  intelligence  about  R7 s  capabilities.  The  approach 
taken  in  this  paper  is  quite  different  as  we  focus  on  developing  deterministic  models  to  address  B’s 
resource  allocation  problems.  The  deterministic  approach  is  justified  in  settings  where  (1)  the  CM 
development  efforts  do  not  involve  a  significant  research  element  and  are  mainly  composed  of  a 
sequence  of  engineering  stages  whose  durations  can  be  forecasted  with  reasonable  accuracy  and 
(2)  when  there  are  reliable  intelligence  reports  regarding  R' s  capabilities,  intentions  and  possible 
hostile  actions. 

The  main  contribution  of  this  paper  is  in  extending  the  operational  situation  described  in  [8]  in 
three  ways:  (a)  assuming  arbitrary  CM  development  policies  (not  necessarily  parallel  or  sequential); 
(b)  introducing  temporal  budget  constraints,  which  are  quite  realistic  in  defense  contracting;  (c) 
allowing  for  a  wide  variety  of  “inconsistent”  CMs  in  the  sense  that  a  certain  CM  may  be  more 
effective  against  weapon  I  than  weapon  II,  while  the  reverse  is  true  for  another  CM.  Also,  unlike 
the  continuous  investment  levels  considered  in  [8],  the  formulation  presented  herein  restricts  the 
investment  levels  to  a  finite  number  of  discrete  values.  Similar  discretization  was  implemented  by 
[2]  to  analyze  investment  levels  among  alternative  projects  related  to  the  natural  gas  industry  in 
the  US. 

We  model  the  decision  problem  of  B  as  a  variant  of  a  resource-constrained  shortest  path  (RCSP) 
problem,  where  the  constraints  capture  global  or  temporal  budgetary  constraints.  The  RCSP 
problem  is  known  to  be  NP-complete,  in  the  ordinary  sense  -  see  [6],  and  we  show  that  our  variant 
is  also  NP-complete.  RCSP  problems  have  been  addressed  by  many  authors  including  [1],  [4], 
[10]  and  [13].  In  particular,  RCSP  problems  of  limited  size  can  be  solved  through  special-purpose 
algorithms  such  as  those  developed  in  [3]  and  [7],  or  through  efficient  general-purpose  algorithms 
available  in  commercial  optimization  software  packages. 


2 


To  demonstrate  the  potential  usefulness  of  our  RCSP  models  and  analyze  their  robustness  to 
small  data  perturbations  we  conducted  an  extensive  computational  study  in  which  we  employed 
the  solver  in  the  MOSEK  optimization  package.  This  solver  was  proven  to  be  quite  efficient  for 
realistically  sized  problem  instances  of  our  RCSP  models. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2  we  introduce  notation  and  state 
the  decision  problem  formally.  In  Section  3  and  4  we  formulate  several  variants  of  the  problem, 
addressing  single  and  multiple  weapons  and  CMs,  as  constrained  network  optimization  models. 
Section  5  demonstrates  the  usefulness  of  the  models  by  presenting  the  results  of  extensive  numerical 
experiments.  Finally,  Section  6  suggests  some  directions  for  future  research. 

2  Problem  Formulation 

The  weapons  that  R  develops  are  indexed  byte  £  W  =  {1,...,|W|}.  For  each  w,  let  sw  denote  the 
time  when  weapon  w  becomes  operational;  the  sw’s  are  obtained  or  estimated  by  B’s  intelligence 
agencies.  In  particular,  sw  =  0  means  that  weapon  w  is  already  operational  at  time  0.  If  -sw  >  0, 
then  weapon  w  does  not  contribute  to  the  damage  inflicted  on  B  until  time  sw.  We  consider  a 
finite  time  horizon  of  length  T.  Without  loss  of  generality,  we  assume  that  T  is  large  enough  such 
that  all  weapons  would  become  operational  before  time  T,  and  that  the  weapons  are  indexed  in 
increasing  order  of  the  sw’s,  that  is,  0  <  si  <  ...  <  S|w|  <  T.  Absent  any  CM,  the  damage  rate 
per  unit-time  inflicted  by  weapon  w  on  B  is  d$  >  0.  Absent  any  weapon,  the  damage  rate  per 
unit-time  is  0,  independent  of  the  available  CMs. 

To  mitigate  the  effect  of  R' s  weapons,  B  develops  CMs.  These  CMs  are  indexed  by  m  £  M.  = 
{1, . . .  ,  |A!|}  and  we  use  the  notation  CM,n  to  refer  to  the  mth  CM.  For  each  w  and  rn,  let  d™  >  0 
be  the  damage  rate  caused  by  weapon  w  when  only  CMm  is  operational.  When  a  set  of  CMs  is 
available,  their  effect  is  not  cumulative  -  the  damage  rate  of  weapon  w  is  determined  by  the  most 
effective  CM  that  is  available  at  that  time.  While  in  some  cases  there  may  be  cumulative  effects 
of  CMs,  e.g.,  when  one  CM  is  a  detection  device  and  another  CM  is  a  neutralization  device,  we 
focus  in  this  paper  on  a  single  family  of  CMs  (e.g.,  interception  systems  or  bomb  neutralization 
systems)  that  evolves  and  improves  over  time  and  whose  members  differ  in  their  capabilities.  So, 
when  a  set  0  ^  M  C  Jv[  of  CMs  is  available,  the  damage  rate  by  weapon  w ,  if  operational,  is 
d m  =  minm€M  d™ ;  when  no  CM  is  available,  the  damage  rate  is  d^  =  d^.  We  will  find  it  useful  to 
apply  the  notation  dM  =  {d\f,  •  •  • ,  d^)  £  for  M  CM.  drn  =  d{m}  for  m  £  M  and  do  =  d$. 
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The  damage  rate  caused  by  a  group  of  weapons  is  represented  by  a  monotonically  increasing 
function  D  :  RIwI  — >  R,  which  converts  damage  rates  of  individual  weapons  into  a  total  damage 
rate.  By  monotonically  increasing  function  D(x i, . . .  >xiw\)  we  mean  that  if  xw  >  yw  for  all  w  = 
1, . . . ,  \W\  then  D(xi, . . .  ,x\w\)  >  D(yi, . . .  ,y\w\ )  and  if  xw  >  yw  for  all  w  =  1, . . . ,  \W\  then 
D(x  i,...,X|iyi)  >  D(yi, . , , ,  y\w\).  For  example,  D(-)  can  be  the  sum  of  the  individual  damage 
rates  or  their  maximum.  When  a  weapon  w  £  W  is  not  operational,  its  individual  contribution 
to  the  total  damage  is  set  to  0  when  evaluating  the  function  D(-).  Given  a  set  IT  C  W,  let 
Iw  G  {0, 1}IWI  be  the  indicator  vector  of  W,  i.e.,  (Iw)w  =  1  if  w  £  W  and  (IW)w  =  0  otherwise. 
The  damage  rate  caused  by  the  set  W  of  operational  weapons  when  the  set  of  available  CMs  is  M 
is  then  expressed  by 

dZ  =  D(IW  o  dM), 

where  “o”  stands  for  the  Hadamard  product.  That  is,  for  vectors  x,y  G  the  Hadamard 

product  is  xo  y  £  Rlwl  with  (io  y)w  =  xwyw.  In  particular,  as  each  d “  represents  the  damage  rate 
of  a  certain  weapon  w  in  the  presence  of  CMm,  we  have  that  D[I^W^  o  dm]  =  dl^.  Also,  D( 0)  =  0. 

It  is  assumed  that  each  CM  can  be  developed  at  any  one  of  several  levels  of  intensity  that  are 
indexed  by&:G/C  =  {l,...,|/C|}.2  A  higher  intensity  level  of  development  has  two  effects:  first,  the 
development  time  of  the  CM  is  shorter  and  therefore  it  becomes  operational  sooner,  and  second, 
the  associated  cost  is  higher.  We  assume  that  the  intensity  of  developing  a  CM  does  not  affect  its 
effectiveness.  Intensity  0  indicates  no  development.  For  m  €  M.  and  k  G  /C,  let  t1^  >  0  denote 
the  time  it  takes  to  complete  CMm  when  developed  at  intensity  level  k  and  let  cj,  >  0  be  the 
corresponding  cost.  As  higher  intensity  levels  are  associated  with  shorter  development  times  and 
higher  costs,  we  index  the  intensities  so  that 

(1)  ^  iml_1  <  -  <  tin  and  cm'  >  4 £l_1  >  -  >  clm  for  each  m  €  M. 

To  avoid  degenerate  situations,  we  assume  throughout  that  there  are  no  ties  among  the  d™ ’ s,  t^’s 
and  c^j’s. 

The  problem  that  B  faces  is  to  decide  which  CMs  to  develop,  at  what  times  to  start  development 
and  at  what  intensity  levels.  The  goal  is  to  minimize  the  cumulative  damage  over  the  time  horizon 
subject  to  budgetary  constraints.  The  simplest  budgetary  constraint  is  a  global  one  where  an  upper 
bound,  say  C.  is  prescribed  on  the  total  funds  that  can  be  spent  in  developing  the  CMs.  In  this 
case,  the  entire  budget  is  available  at  t  =  0.  Since  there  is  no  reason  to  defer  the  development  of 
2The  index  sets  of  the  intensities  of  the  different  CMs  are  the  same  only  for  notational  convenience. 
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any  CM,  the  starting  time  for  all  the  CMs  occurs  at  time  t  =  0.  However,  budgetary  constraints 
may  also  be  temporal  where  there  is  an  upper  bound  on  expenditure  in  a  certain  time  period.  In 
such  cases,  starting  times  of  the  development  of  the  CMs  become  decision  variables. 

3  Consistent  CMs 

In  this  section  we  analyze  variants  of  the  resource  allocation  problem  described  in  Section  2  under 
a  consistency  assumption,  formally  introduced  in  Subsection  3.2.  The  problems  are  formulated  as 
constrained  network  optimization  models,  which  facilitate  their  efficient  solution.  In  Subsection 
3.1  we  consider  the  case  of  a  single  weapon  (which  is  trivially  consistent)  under  a  global  budgetary 
constraint.  In  Subsection  3.2  we  extend  the  analysis  to  the  case  of  multiple  consistent  weapons  and 
in  Subsection  3.3  we  address  temporal  budgetary  constraints. 

3.1  Single  weapon 

Here  we  assume  that  there  is  a  single  weapon  and  no  temporal  budget  constraints;  in  this  case,  the 
convention  that  dff,  =  D[I^ odm\  for  m  G  M.  assures  that  D(-)  is  the  identity,  thus,  d]M  =  D(dM )  = 
dM  for  every  M  C  Jv[ .  Henceforth,  in  this  sub-section  we  suppress  the  index  w  =  1.  Also,  we 
rank  the  CMs  by  their  effectiveness  (against  the  single  weapon),  that  is,  do  >  d\  >  . . .  >  d\j^\  >  0. 
As  there  are  no  temporal  budget  constraints,  it  can  be  assumed  that  the  development  of  all  CMs 
starts  at  0.  We  first  assume  that  the  single  weapon  is  already  operational  at  time  0. 

The  following  definitions  and  notation  are  used  throughout  the  paper.  A  CM- development- 
policy,  henceforth  called  simply  a  policy ,  is  a  set  of  CMs  along  with  corresponding  develop¬ 
ment  intensities.  Thus,  a  policy  7r  is  represented  by  a  set  of  pairs  (rn,  k )  where  the  values 
of  m  are  distinct.  We  order  the  pairs  in  it  by  their  first  coordinate  such  that  the  sequence 
(mi,  k\),  (m2,  kf),  •  •  • ,  (mp,  kp)  satisfies  m\  <  m2  <  ...  <  rnp.  Clearly,  there  is  no  point  in  de¬ 
veloping  a  CM  that  becomes  available  after  a  more  effective  CM  is  already  operational  because 
such  a  CM  cannot  reduce  the  damage  caused  by  the  weapon.  Consequently,  one  can  restrict  atten¬ 
tion  only  to  policies  n  that  have  the  following  property:  if  (m,  k)  and  (mr ,  k')  with  rn  <  m'  are  in 
7 r,  then  <  t^, .  We  refer  to  policies  that  satisfy  this  condition  as  effective  policies. 

A  path  in  a  graph  is  an  ordered  set  of  vertices  where  each  consecutive  pair  is  an  edge.  A  path 
that  starts  at  vertex  a  and  ends  at  vertex  b  is  referred  to  as  an  (a-b)-path.  For  a  path  a,  V(a)  is 
the  set  of  vertices  in  a,  excluding  the  end  vertices,  and  E(a)  is  the  set  of  the  corresponding  edges. 
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Figure  1:  A  FCDS  graph  for  a  single  weapon 

Effective  policies  can  be  represented  by  paths  in  the  acyclic  (directed)  graph  G'  =  ( V',Er ), 
where 

(2)  V'  =  {(m,  k)  :  m  G  M  and  k  £  1C} 
and 

(3)  E'  =  {((m,  k),  (■ m k '))  :  m  <  m!  and  tkm  <  t^,}. 

Next,  we  augment  the  graph  G'  with  a  single  origin  and  a  single  destination,  denoted  O  and  D, 
respectively,  and  with  edges  (O,  D ),  (O,  (to.  A;))  and  ((m,  fc),  D)  for  every  (m,  fc)  G  Eh  The  resulting 
augmented  graph  is  referred  to  as  the  Feasible  CM  Development  Schedule  (FCDS)  graph  and  its 
vertex  and  edge  sets  are  denoted  V  and  E,  respectively.  It  is  convenient  to  represent  the  vertices 
of  a  FCDS  graph  on  the  interval  [0,  T]  with  vertices  of  V'  represented  by  their  corresponding 
completion  times  t^,  O  corresponding  to  0  and  D  corresponding  to  T.  In  this  representation,  all 
edges  have  orientation  from  left  to  right,  assuring  that  the  FCDS  graph  is  acyclic  and  therefore  its 
paths  are  simple.  Figure  1  presents  a  situation  where  developing  CM2  at  intensity  2  takes  less  time 
than  developing  CMi  at  intensity  1  and  therefore  ((1, 1),  (2,  2))  is  not  an  edge  in  the  corresponding 
FCDS  graph.  But,  if  the  same  intensity  level  is  applied  to  CMi  and  CM2,  then  CMi  is  completed 
before  CM2  -  hence,  ((1, 1),  (2, 1))  and  ((1,  2),  (2,  2))  are  edges  in  the  FCDS  graph. 

The  discussion  above  implies  that  we  have  a  one-to-one  correspondence  between  the  set  of 
effective  policies  and  the  set  of  (O-D)-paths  in  the  FCDS  graph.  For  example,  in  Figure  1,  the 
effective  policy  {(1,  2),  (2, 1)}  is  represented  by  the  path  ((O,  (1,  2)),  ((1,  2),  (2, 1)),  ((2, 1),D))  but 
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there  is  no  path  for  the  non-effective  policy  {(1, 1),  (2,2)}.  We  define  two  transformations  denoted 
policy(-)  and  path(-),  which  convert  paths  to  effective  policies  and  vice  versa.  These  transformations 
are  inverses  of  each  other,  i.e.,  if  cr  is  a  path  in  a  FCDS  and  n  is  a  policy,  then  policy  (a)  =n  if 
and  only  if  path(n)  =  cr;  in  particular,  policy  [path(ir)\=TT  and  path  [policy(a)]=a .  Although  the 
definitions  of  policy(-)  and  path(-)  may  seem  trivial,  they  will  be  particularly  useful  in  Section  4. 

The  FCDS  graph  has  |A4||/C|  +  2  vertices  and  at  most  |/C|2(^)  +  2|A4||/C|  +  1  edges.  The 
bound  on  the  number  of  edges  is  attained  when  the  (time)  intervals  {[4n ,  t]m]  :  m  €  A4}  are 
pairwise  disjoint;  that  is,  when  the  development  times  of  the  various  CMs  are  highly  variable.  The 
actual  number  of  edges  of  the  FCDS  graph  depends  on  the  values  and,  in  general,  the  graph 
may  be  sparse  because  of  overlaps  of  the  aforementioned  time- intervals.  That  is,  the  completion 
time  of  some  highly  effective  CMs  developed  at  high  intensity  levels  may  be  shorter  than  that  of 
less  effective  CMs  that  are  developed  at  a  low  intensity  level. 

A  nonempty  effective  policy  ir,  represented  by  a  path  ( O ,  (mi,  k{),  (m2,  fc2),  •  •  • ,  (mp_i, 
kp- 1),  (rrip,  kp),  D)  of  the  FCDS  graph,  defines  a  partition  of  [0,  T]  into  the  time  intervals  [0,  tfy  J,  [i*}  ,  t%*], . . ., 
[tmp-!,tmp\i  [tmp,T].  During  the  time  interval  [tm, , ]  the  available  protection  against  the 
weapon  is  that  of  CMm.,  that  is,  the  damage  rate  inflicted  by  the  weapon  is  dmj.  The  total 
damage  during  that  period  is  dmj  (t-mf+i  —  tmj )  •  Similarly,  the  total  damage  inflicted  during  the 
interval  [0,  ]  is  d^t^  and  the  total  damage  inflicted  during  the  interval  [tmp,  T]  is  dmp{T  —  tmp). 

Thus,  by  assigning  to  each  edge  e  of  the  FCDS  graph  a  damage  value  de  given  by 

~  Am)  if  e  =  ((m,  k),  (m',  k')) 
do(t; -  0)  if  e  =  (O,  (m',  k')) 

(4)  de  = 

dm(T  -  t^)  if  e  =  ((m,  k),  D) 
do(T-  0)  if  e=(0,D), 

the  total  damage  inflicted  in  the  time  interval  [0,  T]  when  an  effective  policy  7 r  is  implemented  is 
expressed  by 

(5)  c?(tt)  =  de- 

e^.E[path(n)] 

The  next  lemma  records  an  opposite  triangular  inequality  that  the  de' s  (the  damage  values) 
satisfy,  which  we  will  use  later  in  Section  4. 

Lemma  1  Suppose  ((mi,  k\),  (m2,  k2)),  ((m2,  k2),  (m3,  /c3))  €  E' .  Then  ((mi,  k\),  (m3,  fe3))  €  E' 
and 

(6)  ^((mi,fci),(m3,fc3))  >  d((rn1,ki),(m,2,k2))  E  4(m2A2)tm3,b))’  ^ 
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Proof:  Trivially,  the  assumptions  imply  m±  <  m .2  <  m3  and  t^f  <  t^2  <  t^3.  To  verify  (6),  note 

that  d((milfc1)l(m3,fc3))  =  dmi  (^m3  —  =  dmi  (tm.3  ~  ^m2)  +  ^mi  (^m2  ~~  ^roj  >  ^2(^1113  —  ^m2)  ~b 

dmi(t%2  =  d((miM)t(m2M))  +d((m2j fc2),(m3,fc3))  (the  last  inequality  follows  from  dmi  >  dm2). 

□ 


The  cost  ce  of  an  edge  in  the  FCDS  graph  represents  the  resources  needed  for  developing  the 
CM  corresponding  to  the  end-vertex  of  that  edge,  that  is, 


(7) 


Cp  — 


if  the  end- vertex  of  e  is  (m,  fc) 


if  the  end- vertex  of  e  is  D. 


The  cost  of  implementing  an  effective  policy  n  is  then  expressed  by 


(8)  c(vr)  =  Y  Ce~ 

e£E[path(  7r)] 

Due  to  (7),  (8)  and  the  one-to-one  correspondence  of  effective  policies  and  paths  in  the  FCDS 
graph,  the  problem  of  selecting  the  effective  policy  that  minimizes  the  damage  subject  to  budget 
constraint  C  reduces  to  the  problem  of  finding  an  (O-D)-path  in  the  FCDS  graph  that  minimizes 
the  cl- length  (5)  subject  to  the  c-length  (8)  being  bounded  by  (7,  that  is, 

(9)  min  )  de 

(7  ‘ 

e£E(a) 

S.t.  ce  <  C 

eEcr 

o  is  an  (O-D)-path  in  the  FCDS  graph. 


Formally, 

Proposition  1  1 n  the  case  of  a  single  weapon,  the  vertices  ( m ,  k)  on  a  path  that  is  optimal  for  (9) 
determine  an  (effective)  optimal  policy  that  minimizes  the  damage  that  R  causes  B  subject  to  the 
budget  constraint. 

Unfortunately,  the  next  result  shows  that  (9)  is  theoretically  hard. 

Lemma  2  The  constrained  shortest  path  problem  on  a  FCDS  graph  is  NP  complete. 

Proof:  See  the  Appendix. 

We  next  relax  the  assumption  that  the  (single)  weapon  is  available  at  time  0  and  assume  that 
it  becomes  operational  at  time  0  <  s  <  T.  In  this  case,  the  total  damage  inflicted  during  a  time 


interval  [a,  b]  in  which  only  CMm  is  available  is  given  by  dm[max{6,  s}  —  max{a,  s}],  and  the  same 
formula  applies  with  m  =  0  if  no  CM  is  available.  It  follows  that  the  total  damage  inflicted  by  the 
weapon  when  an  effective  policy  n  is  used  is  the  length  of  the  corresponding  path  where  edge-lengths 
are  given  by  the  modification  of  (4)  obtained  by  replacing: 

— >■  maxjt*,,  s}  for  each  (m,  k)  G  V 
0  -*•  s. 

No  change  is  needed  in  the  representation  of  the  cost  associated  with  developing  effective  poli¬ 
cies.  It  follows  that  the  problem  still  reduces  to  the  constrained  shortest  path  problem  of  (9)  and 
Proposition  1  extends  to  the  case  where  s  /  0. 

3.2  Multiple  weapons  and  consistent  CMs 

Consider  the  case  where  there  are  multiple  weapons  threatening  B.  We  say  that  the  CMs  are 
consistent  if  the  rankings  of  their  effectiveness  against  all  the  weapons  coincide.  Formally,  this 
means  that  the  CMs  can  be  indexed  so  that 

(11)  d,Q  >  di  >  ...  >  d^M  |  for  each  w  G  W. 

In  particular,  for  any  subset  M  C  At,  argmin m&M  is  invariant  of  w,  implying  that  for  every 

W  C  W 

(12)  djtf  =  D(IW  o  cIm)  =  min  D(IW  o  dm). 

m£M 

The  consistency  property  is  applicable  when  the  weapons  are  similar  (e.g.,  various  types  of  roadside 
IEDs)  and  the  differences  among  the  CMs  are  only  manifested  in  the  extent  of  their  damage 
reduction.  Henceforth,  in  this  section,  we  assume  that  the  CMs  are  consistent.  Also,  no  temporal 
budget  constraints  are  imposed  and  the  development  of  all  CMs  starts  at  0. 

Recall  that  a  policy  is  a  collection  of  pairs  (m,  k)  with  distinct  values  of  m.  The  definition  of 
effective  policies  given  in  Subsection  3.1  relies  on  the  ranking  of  the  effectiveness  of  the  CMs.  Since 
the  rankings  with  respect  to  all  weapons  are  the  same,  it  follows  that  if  CMm  is  more  effective 
than  CM„,/  with  respect  to  one  weapon  then  this  effectiveness  dominance  applies  to  all  weapons. 
Consequently,  if  a  policy  is  effective  with  respect  to  one  weapon,  then  it  is  effective  with  respect  to 
all  weapons. 

Again,  we  start  off  with  the  assumption  that  all  the  weapons  are  operational  at  time  0.  The 
damage  inflicted  by  the  weapons  when  an  effective  policy  n  is  implemented  is  a  monotone  increasing 
function  of  the  damages  inflicted  by  the  individual  weapons.  Specifically,  and  similarly  to  the 
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analysis  in  Subsection  3.1  (see  (4)),  the  total  damage  inflicted  during  the  time  interval  [0,  T]  when 
an  effective  policy  ir  is  implemented  is  expressed  by  (5),  where 

D{IW  o  dm )(t%l,  -  J  if  e  =  ((m,  k),  (m1,  k')) 

D(IW  o  d0)(C  -  0)  if  e  =  (O,  (m* ,  k')) 

D(IW  o  dm){T  -  if  e  =  ((m,k),D) 

D(Iwod0)T  if  e  =  (0,D). 

The  expression  for  the  cost  c(n)  of  implementing  an  effective  policy  n  remains  unchanged  and  is 
expressed  by  (7)-(8).  As  in  Subsection  3.1,  it  follows  that  the  problem  of  selecting  the  best  effective 
policy  reduces  to  the  constrained  shortest  path  problem  of  (9).  Formally, 

Proposition  2  In  the  case  of  multiple  weapons  and  consistent  CMs,  the  vertices  (m,  k )  on  a  path 
that  is  optimal  for  (9)  (with  the  de  values  given  by  (13))  determine  an  optimal  effective  policy  that 
minimizes  the  damage  caused  by  R  to  B  subject  to  the  total  budget  constraint. 

Next  assume  that  some  weapons  are  not  operational  at  time  0  and  weapon  w  becomes  available 
at  time  sw  >  0.  Recall  (from  Section  2)  that  the  weapons  are  indexed  in  increasing  order  of  the 
Sw  s,  i.e. ,  0  <  si  <  ...  <  S|w|  <  T;  thus  potential  sets  of  operational  weapons  that  B  may  encounter 
are  {1,  where  w  G  W.  Using  (12),  we  next  observe  that  B  can  continue  to  restrict  attention 

to  effective  policies.  We  adopt  the  convention  that  [a,  b]  =  0  if  b  <  a  and  r+  =  max{r,  0}  for  a  real 
number  r.  Now,  if  an  effective  policy  n  is  implemented  and  ((m,  k),  (m7,  k'))  €  E\path(n)\,  then  for 
m  G  M.  and  w  G  W,  damage  rate  dm’  is  inflicted  during  the  time  interval  t^,]  (~l  [sw,  sw+i]  = 
[minj^/jS^+il^axl^jSTO}],  whose  length  is  [min{t^,,  sw+i}  -  max{t^,  sw}]+;  when  the  edge 
emanates  from  O,  if),  is  replaced  by  0  and  when  it  terminates  at  H,  t^,  is  replaced  by  T.  So,  for 
e  G  E  and  w  G  W,  let 

[min{f^,,  s^+i}  -  max{f^,  sw}]+  if  e  =  ((m,  k),  {m! ,  k')) 

[min{t^;,  sw+i}  -  sw]+  if  e  =  (O,  (m',  k')) 

[min{T,  sw+i}  -  maxjt^,  sw}]+  if  e  =  ((m,  k),D) 

[min{T,  sw+i}  -  sTO]+  if  e  =  (O,  D). 

°  dm]re \w  if  e  emanates  from  (rn,  k) 

12weW  o  do[Te\w  if  e  emanates  from  O. 

The  total  damage  inflicted  when  an  effective  policy  it  is  used  is  then  the  length  of  the  corresponding 
path  where  edge-lengths  are  given  by  (15).  No  changes  are  needed  in  the  representation  of  the  cost 


(14)  re\w  =  < 


and 

(15)  de  = 


(13)  de  =  { 
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associated  with  the  effective  policy  n.  So,  the  problem  is  still  reduced  to  the  constrained  shortest 
path  problem  of  (9)  and  Proposition  2  extends  to  the  case  where  the  sm’s  are  not  necessarily  0. 

3.3  Multiple  weapons  and  consistent  CMs  with  temporal  budget  constraints 

In  this  subsection  we  consider  the  situation  examined  in  Subsection  3.2  with  additional  constraints 
that  restrict  periodical  expenditures.  We  specify  H  time-intervals  (subsets  of  [0,  T])  as:  I\  = 
[T,  ,  T\\.  I2  =  [T_2,  T2],  ■  ■  ■ ,  Ih  =  [Hhi  T h\,  where  T_h  <Th  and  assume  that  there  is  a  bound  Ch 
on  the  expenditure  during  time  interval  Ih,  for  each  h  =  1,2 in  addition  to  the  global 
budget  constraint.  The  time  intervals  are  not  necessarily  disjoint.  There  are  two  special  cases  of 
particular  interest.  In  the  first  case,  the  intervals  Ih  are  disjoint  and  thus  partition  [0, T\.  In  this 
case  the  periodic  budget  constraints  represent  strict  cash  flow  constraints  (e.g.,  typical  to  the  US 
government  budgeting  rules)  in  which  excess  funds  in  one  period  cannot  be  utilized  in  the  next 
period.  Clearly,  for  the  global  constraint  to  be  meaningful,  the  data  must  satisfy  J2h=i  Ch  >  C. 
The  second  case  involves  relaxed  cash  flow  constraints  where  budget  overflows  are  allowed  to  be 
used  in  future  periods.  In  this  case,  the  time  epochs  0  =  T\  <  T2  <  ...  <  Th  <  Th+ 1  =  T  are 
given  and  Ih  =  [0,  for  h  =  1  Here,  for  the  data  to  be  meaningful  it  must  satisfy 

C±  <  C2  <  ■■■  <  Ch  <  C  (the  constraint  corresponding  to  [0,  Th+i]  is  left  out  as  it  is  the  global 
budget  constraint). 

We  assume  that  costs  are  incurred  continuously  and  at  a  uniform  rate.  Also,  the  development 
of  CMs  is  carried  out  without  planned  interruptions  (an  assumption  that  is  quite  reasonable  since 
in  reality,  disrupting  a  project  may  incur  high  set-up  cost  when  development  is  resumed).  In 
the  presence  of  bounds  on  periodic  expenditures,  the  time  at  which  the  development  of  each  CM 
starts  becomes  a  decision  variable.  To  address  these  additional  decision  variables,  we  extend  the 
definition  of  “intensity”  and  refer  to  “plans”  of  CM  development  projects.  Each  plan  consists  of  a 
pair  ( k,r ),  where  k  is  the  intensity  of  developing  the  CM  and  r  is  the  start  time.  For  simplicity  of 
exposition,  we  focus  on  the  case  where  neither  the  effectiveness  of  a  CM,  nor  the  cost  and  duration 
of  developing  it  are  affected  by  r;  relaxation  of  these  assumptions  is  briefly  discussed  at  the  end  of 
this  subsection. 

As  discussed  above,  when  there  are  no  temporal  budgetary  constraints,  there  is  no  reason  to 
defer  the  development  of  any  CM  and  the  t^n  values  represent  both  the  development  duration  and 
the  completion  time.  But,  this  identity  does  not  hold  when  the  starting  times  of  developing  the 
CMs  are  not  0.  Specifically,  if  CMm  is  developed  using  plan  ( k,r ),  then  its  completion  time  is 
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tm’T'>  =  r  +  tjjj.  To  avoid  irrelevant  situations,  we  assume  that  tm <  T  for  all  m,  k  and  r. 

A  policy  77  is  now  defined  as  a  set  of  triplets  (m,  k,  r)  with  distinct  values  of  m  and  effective 
policies  are  defined  in  terms  of  the  tm'T ^  (the  completion  times)  rather  than  the  t^’s  (the  duration 
times) .  Consider  the  modification  of  the  FCDS  graph  where  vertices  and  edges  are  defined  in  terms 
of  the  triplets  (m,  k,  r)  and  the  completion  times  tm’T ^  instead  of  the  pairs  (rri.  k )  and  the  duration 
times  t% j.  We  refer  to  the  resulting  graph  as  the  Timing- Feasible  CMs  Development  Schedule  (T- 
FCDS)  graph.  Note  that  if  the  development  of  all  CMs  can  start  at  any  one  of  q  potential  time 
periods  and  /C  is  the  set  of  potential  intensity  levels,  then  the  T-FCDS  graph  has  g|A4||/C|  +  2 
vertices  and  at  most  (q\IC\)2  (^)  +  2q,|A/1||/C|  +  1  edges.  However,  unlike  the  FCDS  graph,  the 
T-FCDS  graph  may  be  quite  dense  in  real-world  applications.  Since  an  effective  policy  is  such 
that  for  any  two  nodes  ( m,k,T )  and  ( m',k',T '),  m  <  m',  in  the  T-FCDS  graph  we  have  that 
tm’  '  <  tm, ’  ,  it  follows  that,  as  before,  there  is  a  one-to-one  correspondence  between  the  set  of 

effective  policies  and  the  set  of  (O-D)-paths  in  the  T-FCDS  graph. 

Suppose  that  the  development  plan  of  CMm  is  ( k,r ).  During  the  development  period  [r,  tm'T'>] 
a  cost  frm  =  -njr  per  unit  time  is  incurred.  The  total  expenditure  on  CMm  during  the  time  interval 
Ih  =  [Th,Th]  is  then  the  length  of  the  (possibly  empty)  interval  J/jflfr,  tm’T'>]  times  the  per  unit-time 
cost  c*,.  We  note  that  Ih  fl  [r,  tm  =  [max{Tft,  r },  min{T/j,  tm’T^  }].  Consequently,  the  expenditure 
during  the  time  interval  Ih  on  CMm  is  c^jfminjT/j,  tm’^}  ~  nrax{Tft,  r}]  +  .  For  each  edge  e  of  the 
T-FCDS  graph  and  h  =  0, 1, . . . ,  H,  let 


(16) 


c^[mm{Th,tmT)  }-  max{T^,  r}]+  if  e  terminates  at  (■ m,k,r ) 
0  if  e  terminates  at  D. 


The  total  cost  associated  with  effective  policy  n  during  the  time-interval  Ih  is  then  expressed  by 
^2eeE]path(ir)]  ce  and  this  sum  is  subject  to  the  corresponding  temporal  budget  constraint  Ch. 

We  next  consider  the  total  cost  associated  with  implementing  effective  policy  it.  Similarly  to 
the  case  of  a  single  weapon  presented  in  Subsection  3.1,  here  the  total  cost  is  expressed  by  (8), 
with  the  ce’s  given  by  (7),  except  that  ( m,k )  is  replaced  by  (m,  k,r). 

Now,  consider  the  total  damage  associated  with  an  effective  policy  77.  First,  assume  that  all 
weapons  are  operational  at  time  0.  For  each  edge  e  of  the  T-FCDS  graph,  let  de  be  given  by  the 
variant  of  (13)  in  which  k  and  k1  are  replaced  by  (fe,r)  and  (k',Tr),  respectively.  The  total  damage 
associated  with  policy  77  is  then  expressed  by  (5).  When  weapon  w  becomes  operational  at  time 
sw  >  0  ,  let  de  be  given  by  (14)  and  (15),  with  (fe,r)  replacing  k. 

The  above  discussion  demonstrates  that  the  problem  of  selecting  a  policy  that  minimizes  the 
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total  damage  subject  to  total  and  temporal  budget  constraints,  reduces  to  the  problem  of  finding  an 
(O-D)- path  in  the  T-FCDS  graph  that  minimizes  the  d- length  subject  to  corresponding  constraints 
on  the  c- length  and  the  c,!- lengths,  that  is, 

(17)  min  ^  de 

e(zE[cr\ 

S.t.  ^2  ce  <  C 

eGE[a] 

J2  ce  —  Ch  for  h  =  1, . . . ,  H 

eGE[a] 

a  is  an  (O-H)-path  in  the  T-FCDS  graph. 

Formally, 

Proposition  3  In  the  case  of  multiple  weapons  and  consistent  CMs  with  temporal  budget  con¬ 
straints,  the  vertices  ( m,k,r )  on  a  path  that  is  optimal  for  (17)  determine  an  optimal  effective 
policy  that  minimizes  the  damage  caused  by  R  to  B  subject  to  the  total  and  temporal  budget  con¬ 
straints. 

We  note  that  our  model  can  be  easily  modified  to  capture  situations  where  costs  and  devel¬ 
opment  times  of  CMs  depend  on  the  starting  times  of  their  development;  all  that  is  needed  is  to 
replace  the  parameters  c(:n  and  tff  by  start-time-dependent  counterparts  c^(r)  and  t^(r),  respec¬ 
tively  (in  which  case  fm,T  =  r  +  t^fr)).  Also,  allowing  effectiveness  of  the  CMs  to  depend  on  the 
starting  times  can  be  captured  by  replacing  rn  with  (m,r)  in  the  modification  of  (13). 

4  Inconsistent  CMs 

In  this  section  we  relax  the  consistency  assumption,  allowing  for  one  CM  to  be  more  effective  than 
another  with  respect  to  weapon  w  while  the  reverse  holds  for  weapon  w'  w.  Inconsistency  is 
present,  for  example,  when  the  weapons  of  R  are  not  technologically  or  operationally  similar.  We 
analyze  the  inconsistent  case,  only  when  the  function  Df)  is  the  summation  function,  in  which 
case  D(IXX  o  d\i)  =  Ylw  gw  ^m-  We  start  our  analysis  under  the  assumptions  that  all  weapons  are 
operational  at  time  0  and  there  are  no  temporal  budgetary  constraints. 

Without  the  consistency  assumption,  the  CMs  can  no  longer  be  ranked  uniformly  according  to 
their  effectiveness  against  the  weapons.  Still,  each  weapon  has  its  own  total  order  regarding  the 
effectiveness  of  the  CMs,  which  we  call  re-domination  and  denote  by  -<w ;  thus,  we  write  m  -<w  m' 

if  C'<C- 
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4.1  Policies  and  their  effective  parts 


We  return  to  the  definition  in  which  a  policy  n  is  a  set  of  pairs  (m.  k )  with  distinct  values  of  m. 
A  policy  is  w- effective  if  it  is  effective  in  the  sense  of  the  definition  in  Subsection  3.1  for  weapon 
w.  Given  a  policy  n  and  a  weapon  w  G  W,  the  w-effective  part  of  it,  denoted  irw,  is  the  subset  of 
7T  obtained  by  removing  all  pairs  ( m k')  G  n  for  which  there  exists  (rn.  k )  G  7r  with  m'  -<w  m  and 
tm  <  t-m'-  Since  a  policy  applies  at  most  one  intensity  level  for  each  CM,  the  irffs  must  satisfy  the 
following  “coupling  requirement” : 

(18)  [(m,  k )  G  irw  and  (m,  k')  G  7rw/]  =4>  [k  =  k']. 

The  total  cost  of  implementing  a  policy  is  expressed  by 

(19)  c(tt)  =  cm- 

(m,k)^.TT 

4.2  The  multi-FCDS  graph 

Each  weapon  w  defines  an  FCDS  graph,  denoted  FCDSW,  where  <  in  (3),  with  respect  to  the  CM 
indices,  is  replaced  by  -<w.  The  sets  of  vertices  and  edges  of  the  FCDS^,  graph  are  denoted  Vw  and 
Ew ,  respectively.  Each  Vw  is  a  replica  of  U  =  {0,D}  U  {(m,  k)  :  m  G  M,k  G  /C};  in  particular, 
the  elements  of  Vw  are  denoted  Ow,  Dw  and  ( m,k)w .  Since  the  CMs  are  inconsistent  and  an 
edge  (( m,k)w ,  (m',k')w)  exists  in  Ew  if  and  only  if  m  -<w  m'  and  t ^  <  t^,  ,  it  is  possible  that 
((in,  k)w,  (m',  k')w)  G  Ew  but  ((m,k)w’ ,(m',k')w')  ^  Ew'  for  w'  w  (this  will  happen  if  and 
only  if  f <  t^,,  m  -<w  m'  and  m'  -<w' m).  As  (( m,k)w ,  (m',k')w)  G  Ew  implies  that  m  -<w  m', 
the  FCDSW  graphs  are  acyclic  and  therefore  their  paths  are  simple. 

In  order  to  refer  to  the  FCDS^,  graphs  jointly,  we  define  the  multi-FCDS  graph  (V,  E)  whose 
vertex-  and  edge-sets  are,  respectively,  V  =  U w&y\>Vw  and  E  =  U w£\yEw.  The  FCDSW  subgraphs 
are  disjoint  components  of  the  multi-FCDS  graph,  in  particular,  (Ow-Dw)- paths  of  the  FCDS^, 
graph  are  identified  with  the  (0“'--Dw )-paths  of  the  multi-FCDS  graph.  Figure  2  presents  a  multi- 
FCDS  graph  for  two  weapons.  As  in  Figure  1,  here  too,  t\  <  t\  <  t\  <  t\.  However,  while  CMi  is 
more  effective  than  CM2  with  respect  to  weapon  1,  the  reverse  is  true  for  weapon  2. 

Results  in  Subsection  3.1  show  that  there  is  a  one-to-one  correspondence  of  -w-effective  policies 
and  (O^-D^j-paths  of  the  FCDS^  graph,  which  are  the  (Ow-Dw)-paths  of  the  multi-FCDS  graph. 
Following  the  notation  introduced  in  Subsection  3.1,  the  (Ow-D“)-path  corresponding  to  a  in¬ 
effective  policy  irw  is  denoted  path(irw).  Figure  3  demonstrates  the  tc-effective  parts  of  a  policy. 
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Figure  2:  A  multi-FCDS  graph  for  two  weapons 


Here,  7 r  =  {(1, 1),  (2,  2)};  its  1-effective  part  is  {(1, 1),  (2, 2)}(=  7r),  whereas  its  2-effective  part  is 
only  {(2,2)}. 

For  each  e  G  E,  let 

if  e  =  ((in,  k)w,  (m',  fc')"’) 

if  e  =  (0”,(m',vn 

(20)  de  =  < 

<(T-4)  if  e  =  ((m,  k)w,Dw) 

d%T  if  e  =  (Ow,Dw). 

Similar  to  (5)  in  Subsection  3.1,  we  then  have  that  the  damage  that  weapon  w  €  W  causes  when 
policy  7T  is  implemented  is  expressed  by  d(irw)  =  YleeE\path(irw)]  de-  The  total  damage  associated 
with  the  implementation  of  policy  7r  (caused  by  all  weapons)  is  then  expressed  by 

(21)  d( tt)  =  y  dM  =  Y  Y  de' 

w£\V  uigW  e£E[path(TTw)] 

4.3  Configurations 

Unlike  the  one-to-one  correspondence  of  effective  policies  and  paths  in  Subsections  3. 1-3.3,  the 
relation  between  policies  and  paths  in  the  inconsistent  case  is  more  complex  because  the  7rw’s  of  a 
policy  7r  can  differ  from  each  other.  To  overcome  this  difficulty,  we  define  a  configuration  to  be  a 
collection  a  =  {aw  :  w  €  W},  where  each  aw  is  an  (Ow-Dw)- path  in  the  FCDSW  subgraph  of  the 
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policy  (7t) 

*1 
7^2 

Figure  3:  A  representation  of  the  tc-effective  parts  of  a  policy 
multi-FCDS  graph.  Such  a  configuration  is  called  plausible  if  it  satisfies  the  “coupling  requirement” 

(22)  [(to,  k)w  G  V(aw)  and  ( m,k')w  G  V{aw>)\  =4>  [A;  =  k']. 

Note  that  (18)  implies  that  the  paths  corresponding  to  the  u;-effective  parts  of  a  policy  tt  form  a 
plausible  configuration.  On  the  other  hand,  given  a  plausible  configuration  a  =  {aw  :  w  G  W}, 
we  define  policy(a)  =  'Jw&v\>V{aw).  Finally,  we  point  out  that  policy  tt  consists  of  nodes  from  the 
multi-FCDS  graph  that  may  be  replications  of  (to,  k )  pairs  e.g.,  (m,  k)w  and  (to,  k)w' ,w  =  v/ .  The 
plausibility  of  a  assures  that  policy(a)  is  indeed  a  policy.  The  transformations  n  — >  { path(rrw )  : 
re  G  W}  and  <r  — >■  policy{a)  are  not  inverses  of  each  other.3 
For  a  policy  -it  we  define 

(23)  tt  =  policy ({path(TTw)  :  w  G  W}); 

3For  example,  it  is  possible  to  have  a  policy  n  which  is  not  w  -effective  for  any  w  €  W  -  while  n  determines 
a  configuration,  it  is  not  in  the  image  of  policy(-).  Also,  given  a  configuration  a  and  w  ^  w' ,  it  is  possible  for 
aw  to  contain  a  vertex  (m,  k)  that  is  not  u/-dominated  by  any  vertex  of  awi  -  this  vertex  will  then  appear  in 
path[policy(a)]wi . 


o 
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note  that  jt  C  j r  and 

(24)  c(ir)=  cm>  E  cm  =  c(^); 

(m,fc)G  7r  (m,fc)G7r 

strict  inclusion  7f  C  7t  and  strict  inequality  in  (24)  are  possible  when  7 r  contains  “superfluous” 
(m,  &;)’s.  Also,  for  each  path  aw  of  a  configuration  <7,  ^(cq,,)  C  \policy(cr)\w  and  therefore,  by  the 
opposite  triangular  inequality  (Lemma  1), 

(25)  E  E  d- 

e£E(aw)  e£E{path[policy(cr)w]} 

Combining  (25)  with  (23)  shows  that  for  every  configuration  <7  =  {aw  :  w  G  W}: 

(26)  d\policy(a)\  =  E  E  "«<E  E 

wEW  e£E{path[policy(<j)w]}  W  e£E(crw) 

Figure  4  demonstrates  a  situation  with  strict  inequality  in  (26).  The  figure  illustrates  a  plausible 
configuration  <7  with  a\  =  (O1,  (1,  l)1,  D1)  and  02  =  (O2,  (2,  2)2,  D2)  (the  edges  of  these  paths  are 
marked  in  bold).  The  policy  corresponding  to  a  is  7r  =  {(1, 1),  (2,2)}  and  the  path  corresponding 
to  the  1-effective  part  of  7r  is  (O1,  (2, 2,  )1,  (1,  l)1,  D1),  with  edges  (O1,  (2, 2,  )1)  and  ((2, 2,  )1,  (1,  l)1) 
(marked  by  dotted  arrows)  replacing  edge  (O1,^,  l)1)  of  oi.  Due  to  the  opposite  triangular  in¬ 
equality,  {(1, 1, ),  (2,  2)}  provides  better  protection  against  weapon  1  than  {(1, 1)}  -  the  1-effective 
policy  corresponding  to  a\ . 

4.4  Network  optimization  with  side  constraints 

For  an  edge  e  and  vertex  v  of  the  multiple-FCDS,  we  write  e  <—  v  if  e  emanates  from  v  and  e  — »  v 
if  e  terminates  at  v.  Standard  results  of  network  modeling  imply  a  one-to-one  correspondence 
between  configurations  and  vectors  x  =  (xe)e^E  that  satisfy: 

^2e<—(m,k)w  xe  =  Ee^(m,k)™  xe  f°r  each  m  G  M,k  G  1C  and  W  G  W 

(27)  Ee^Ow  xe  =  1  =  Ee->DW  Xe  for  eaC^  W  €  TV 

xe  G  {0, 1}  for  each  e  G  E. 

In  particular,  if  configuration  a  =  {aw  :  w  G  W}  corresponds  to  x  satisfying  (27),  then 

(28)  E  E  de  =  E deXe' 

w€W  e£E(aw)  eG E 

Given  such  a  solution  x,  let  y{x)  =  [y(a:)m](m,fc)Gt/  be  defined  by 

y(x)m  =  max  %e  for  each  ( m,k )  G  U\ 

w  Z — J 

e—>(m,k)w 
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-  Plausible  configuration 
. . .  Improved  plausible  configuration 

Figure  4:  Demonstrating  opposite  triangular  inequality 
equivalently,  y(x )  is  the  unique  solution  of 

y{x)m  -  xe  for  each  (m>  k)  £  U  and  w  £  W 

y(x)m  <  E^ew  Ee-(m,fc)«  for  each  (m,  fc)  €  C7 

y(x)^  G  {0, 1}  for  each  (m,  k )  €  U. 

We  note  that  the  variable  y(x)%l  gets  the  value  1  if  for  any  weapon  re,  CMm  is  developed  at  level  k 
under  x  when  considering  the  tc-subgraph  of  the  multi-FCDS  graph,  in  all  other  cases  y(x)^n  gets 
the  value  0. 


Consider  a  configuration  a  corresponding  to  the  vector  x  =  ( xe)e&E  that  satisfies  (27).  Evi¬ 
dently,  <7  is  plausible  if  and  only  if  EfcefC  —  1  fo1'  each  m  £  M;  an  equivalent  condition  is 

that  for  some  vector  y  =  (ym){m,k)eU:  x  and  1J  satisfy 


(29) 


Vm  >  E 


e—>(m,k)w 


Vm  —  EidGW  E 


e— >(m,k)w  xz 


E fee/C^/m  —  1 

ykm  €  {0, 1} 


for  each  (rn.  k )  G  U  and  w  £  W 
for  each  (rn.  k )  G  17 
for  each  m  £  M 
for  each  (m,  7)  G  U. 
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Thus,  we  have  a  one-to-one  correspondence  between  the  set  of  plausible  configurations  and  the  set 
of  solutions  of  (27)  and  (29).  Further,  if  (x,y)  satisfying  (27)  and  (29)  corresponds  to  a  plausible 
configuration  a,  then  policy  (a)  =  {( m,k )  :  =  1}  and 

(30)  c[policy(o )]  =  £  ckmykp, 

( m,k)£U 

the  budget  constraint  on  policy{a)  is  then  expressed  by 

(31)  £  ckmykm  <  C. 

(■ m,k)£U 

We  next  show  that  the  problem  of  selecting  an  optimal  policy  reduces  to  solving  the  optimization 
problem: 

(32)  min  )  dexe 

x,  y  J 

eeE 

s.t.  (x,y)  satisfies  (27),  (29)  and  (31). 

Proposition  4  I n  a  setting  that  possibly  involves  inconsistent  weapons,  let  (x*,y*)  be  an  optimal 
solution  of  (32)  and  n*  =  {(m,  k)  G  U  :  (y*)m  =  !}•  Then  tt*  is  a  policy  that  minimizes  the  damage 
caused  by  R  to  B  subject  to  the  total  budget  constraint. 

Proof:  From  the  definition  of  n* , 

(33)  <ko=  E  s  c. 

(■ m,k)eU 

assuring  that  n*  satisfies  the  budget  constraint.  Next,  let  OPT  be  the  optimal  value  of  (32). 
With  a*  =  {af,  :  w  G  W}  as  the  plausible  configuration  corresponding  to  (x*,y*),  we  have  that 
7 t*  =  policy  (a*)]  it  then  follows  from  (26)  and  (28)  that 

(34)  d( tv*)  =  d\policy(cr*)}  <  £  £  =  £  £  dex*e  =  OPT. 

w&W  e£E(a *)  w&V\>  e&E 

To  see  that  7r*  is  optimal  consider  an  arbitrary  policy  n  that  satisfies  the  total  budget  constraint. 
Let  (. x ,  y)  be  the  solution  of  (27)  and  (29)  that  corresponds  to  the  plausible  configuration  o  = 
{pat.h(irw)  :  w  G  W}.  Then  policy(o)  =  tt  and  (30)  and  (24)  imply  that 

£  ckmykm  =  c(tt)  <  c(tt)  <  C, 

(m,k)GU 

assuring  that  (. x,y )  satisfies  (31).  So,  (x,y)  is  feasible  for  (32)  and  therefore  OPT  <  Yhe&E  dexe-  It 
now  follows  from  (34),  (28)  and  (21)  that 

d(TT*)  <  OPT  <  £  dexe  =  £  £  de  =  d( tt).  □ 

e$LE  wEW  e£E[path(7rw)] 
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We  next  consider  the  case  where  weapons  are  not  necessarily  operational  at  time  0  and  for  each 
w  G  W,  weapon  w  becomes  operational  at  time  0  <  sw  <  T .  The  situation  is  handled  by  adjusting 
the  definition  of  the  damage  coefficients  given  in  (20)  (resembling  the  use  of  (10)  to  modify  (4)  in 
Section  3.1).  Specifically,  for  m  G  M,k  G  /C  and  w  G  W,  let  tk{u'  =  max{t^,sw}  and  for  each 
e  G  E,  let 

dZitm'™  ~  if  e  =  ((”»>  k)w,  (m1,  k')w) 

do{tkJw  -  sw)  if  e  =  ( Ow ,  (m',  k')w) 

de  = 

dZ(T-t^w)  if  e  =  {(m,k)w,Dw) 

d$(T-sw)  if  e  =  (Ow,Dw). 

With  this  adjustment,  the  total  damage  inflicted  when  an  effective  policy  ir  is  used  is  the  sum  of 
the  length  of  the  corresponding  paths,  the  decision  problem  of  B  reduces  to  (32)  and  Proposition 

4  applies.  It  is  also  possible  to  generalize  (32)  to  the  case  where  periodic  budget  constraints  are 
imposed  (see  Subsection  3.3).  While  the  model  remains  essentially  the  same,  its  size  increases 
significantly.  The  numerical  experiments  reported  in  the  next  section  apply  to  this  more  general 
case  with  cash  flow  budget  constraints.  That  is,  the  intervals  Ih  partition  [0,  T]  with  time  epochs 
0  =  T\  <  T-2  <  ...  <  Th  <  Th+ i  =  T  and  each  interval  Ih  corresponds  to  the  period  [Tjj,  Th+ 1]. 

5  Numerical  Experiments 

In  order  to  examine  the  applicability  and  robustness  of  our  constrained  network  model,  we  con¬ 
ducted  an  extensive  computational  study.  In  this  study,  we  implemented  a  model  with  inconsistent 
weapons  (see  Section  4),  multi-period  temporal  budget  constraints  (See  Subsection  3.3)  where  the 
temporal  budget  constraints  correspond  to  H  disjoint  intervals  that  partition  [0,  T],  and  D(-)  as 
the  summation  function  (see  the  last  paragraph  of  Subsection  4.4).  A  detailed  formulation  of  the 
model  implemented  in  our  numerical  study  is  given  in  the  Appendix. 

The  case  study  consisted  of  10  base  cases  and  the  generation  and  solution  of  100  instances 
for  each  of  the  base  cases.  Each  instance  represented  a  small  perturbation  of  the  development 
completion  times  of  the  CMs.  The  specific  model  we  address  is  provided  in  the  appendix.  All 
computational  tests  were  carried  out  on  a  laptop  computer  with  a  Genuine  Intel  1  GHz  T2500 
processor  and  1  GB  RAM,  running  the  Red  Hat  Linux  5  operating  system.  The  code  was  built 
with  C++  version  4.1.1  and  linked  with  glibc  2.5.4  and  Mosek  5.0.  Each  one  of  the  10  base  cases 
represented  10  types  of  (inconsistent)  weapons  (|W|  =  10),  10  possible  CMs  (|A4|  =  10),  3  levels 
of  intensity  (|/C|  =  3)  and  4  time  periods  ( H  =  4).  Consequently,  each  integer  programming 
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formulation  has  about  38, 000  variables  and  73, 000  constraints  and  the  corresponding  multi-FCDS 
graph  has  about  120  vertices  (one  vertex  for  each  possible  completion  time  tm’T and  4,000  edges. 
Although  the  problems  are  quite  large,  the  model  proved  to  be  computationally  efficient;  the  average 
running  time  of  each  instance  was  about  15  minutes  (the  shortest  time  was  just  a  few  seconds  and 
the  longest  time  was  a  little  over  4  hours). 

5.1  Parameters  setting 

The  CMs  were  divided  into  3  groups  according  to  their  completion  times:  the  first  group  (consisting 
of  CMs  indexed  by  m  =  1, 2,  3)  had  short  development  times,  the  second  group  (consisting  of  CMs 
indexed  by  m  =  4,  5, 6)  had  medium  development  times  and  the  third  group  (consisting  of  CMs 
indexed  by  m  =  7, 8,  9, 10)  had  long  development  times.  The  time  units  were  measured  in  months 
and  each  time  period  was  set  to  be  12  months.  So,  the  first  period  is  [0, 12],  the  second  [12,  24]  etc., 
thus  T2  =  12,  X3  =  24  and  T4  =  36.  The  time  horizon,  T,  was  defined  to  be  T4  +  maxk,m{tm}- 
For  simplicity,  in  our  computational  testings  we  limit  r,  the  optional  starting  times  of  developing 
the  CMs  to  take  the  values  of  0  =  Ti,T2,T-^  or  T4.  The  completion  times,  the  t^’s,  were  sampled 
uniformly  from  the  data  listed  in  Table  1. 


k 

m 

1 

2 

3 

1,2,3 

4,6,8 

10,12,14 

16,18,20 

4,5,6 

12,15,18 

21,24,27 

30,33,36 

7,8,9,10 

24,28,32,36 

40,  44,  48 

52,56,60 

Table  1:  Values  for  the  t^’s 


For  each  CM,  the  development  completion  times  were  ordered  according  to  the  three  intensities 
such  that  tlm  >  >  t'm-  The  damage  rates,  the  d™’ s,  were  sampled  from  a  trimmed  Normal 

distribution  with  parameters,  fi  =  2,  <7  =  1.5  and  truncated  by  10  and  0  as  upper  and  lower 
bounds,  respectively.  Recall  that  in  the  case  of  inconsistent  CMs,  the  d™  ’ s  do  not  have  the  same 
order  of  the  m’s  for  all  w  €  W;  still,  for  each  w,  dg  is  the  maximum  over  all  d™.  In  addition,  sw, 
the  time  weapon  w  is  expected  to  be  operational,  was  sampled  from  the  Uniform  distribution  on 
the  interval  [0, 15]  and  these  times  were  sorted  in  an  ascending  order.  The  development  costs 
were  sampled  from  the  Uniform  distribution  on  the  interval  [3, 10]  and  for  each  to,  the  c^’s  were 
ordered  to  be  nondecreasing  in  k.  The  budget  temporal  constraints  for  all  time  periods  were  set  to 
be  identical  for  all  h,  so  that  Ch  =  C  for  all  h.  C  was  calculated  as  the  average  costs  multiplied 
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by  twice  the  average  period’s  duration  5,  that  is: 


Ey'  cm 

fce/c  l^meM  0 
\IC\-\  M\  -2'5’ 

V"'//  /  m  ryn  \ 

where  6  =  z^h=l(  J^+1 — — .  The  budget  constraint  for  the  entire  time  period  T,  was  taken  as  C  = 
C'-tf  -0.9. 

5.2  Testing  robustness 


As  indicated  before,  experiment  was  carried  out  in  order  to  test  the  robustness  of  the  solutions  to 
mild  perturbations  in  the  data.  In  particular,  we  investigated  the  sensitivity  to  small  deviations  in 
the  development  completion  times.  First,  the  10  base  cases  were  solved.  We  next  generated  100 
instances  for  each  base  case,  with  each  instance  representing  a  mild  perturbation  of  the  completion 
times  of  the  base  case.  The  instances  were  generated  by  sampling  from  a  trimmed  Normal  distribu¬ 
tion.  Specifically,  for  the  completion  time  of  the  mth  CM  under  the  kth  intensity  we  used  a  Normal 
distribution  (/i,  a)  with  p  =  and  a  determined  in  the  following  way:  for  m  =  1,  2, 3  :  a  =  0.05T4, 
for  m  =  4,  5, 6  :  a  =  0.05  •  24  and  for  m  =  7, 8, 9, 10  :  a  =  0.05  •  44.  The  range  of  the  trimmed 
Normal  distribution  was  [p  —  3er,  p  +  3a].  The  completion  times  were  not  re-ordered  to  be  monotone 
in  the  development  intensity  level.  The  time  intervals  and  the  damage  costs,  which  are  functions 
of  the  completion  times,  were  re-calculated  in  the  same  way  as  in  the  base  case  while  the  Ci  s  were 
not  recalculated.  Each  of  the  100  instances  of  each  base  case  was  solved  separately. 

5.3  Measurements 


We  tested  the  robustness  of  the  solution,  relative  to  each  base  case  j  €  {1, . . .  ,  10}  and  each 
perturbation  i  6  {1, . . . ,  100}.  For  each  j  and  i,  let  PertOptji  be  the  optimal  value  (total  damage) 
of  the  ith  perturbation  of  the  jth  base  problem,  let  ( X°pt ,  Y°pt )  be  an  optimal  solution  of  the  jth  base 
case,  and  let  PertBasisOptji  be  the  total  damage  incurred  if  the  actual  development  completion 
times  are  as  in  the  jith  instance  and  (X°pt,  Y°pt)  is  used.  We  consider  the  following  measure  of 
robustness: 

( PertBasisOptji  —  PertOptji) 

^r‘  PertBasisOptji 

Note  that  (Xjpt,Y°pt)  might  be  an  infeasible  solution  to  the  ith  perturbation  of  the  jth  base  case. 
Infeasibility  might  be  of  two  types  -  violation  of  the  topology  or  violation  of  the  temporal  budget 
constraints.  Topological  violation  can  occur  if  some  of  the  edges  in  the  solution  corresponding  to 
(X°pt\  Y°pt)  do  not  exist  in  the  graph  corresponding  to  the  jith  instance.  Such  violations  were 
handled  by  direct  calculation  of  the  actual  value  of  the  solution  of  the  jith  instance.  We  did  not 
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calculate  the  pji  s  for  instances  that  were  determined  infeasible  due  to  violations  of  the  temporal 
budget  constraints. 


5.4  Numerical  results 

Table  2  and  its  histogram  in  Figure  5  illustrate  the  frequencies  of  the  p  values  obtained  from 
792  feasible  solutions  out  of  the  1000  runs  (in  208  runs  there  were  mild  violations  in  the  budget 
constraints  as  indicated  in  Tables  3  and  4). 


Value 

0 

.004 

.009 

.013 

.018 

.022 

.027 

.031 

.036 

.040 

more 

Frequency 

456 

127 

86 

39 

25 

24 

16 

12 

5 

2 

i 

Table  2:  Frequency  of  the  p  values  in  the  original  runs 


Count 

Average 

STD 

Max 

Min 

792 

0.0040 

0.0073 

0.0446 

0 

Table  3:  p  values  with  the  original  temporal  budget  constraints 


We  further  examined  the  effect  of  the  temporal  budget  violations  in  the  following  way.  For 
each  infeasible  perturbation,  we  re-ran  the  model  with  an  increase  of  5%  in  each  violated  temporal 
budget  constraint.  The  purpose  of  this  exercise  was  twofold  -  first,  to  see  whether  a  relatively 
small  increase  in  the  budget  (5%)  is  sufficient  to  restore  feasibility  and  second,  to  see  what  effect 
will  that  increase  have  on  the  p  values.  Table  4  and  its  histogram  in  Figure  6  provide  the  resultant 
p  values  after  the  5%  increase  was  implemented  where  necessary  and  Table  5  provide  the  new 
summary  statistics.  Comparing  Table  3  to  Table  5  we  observe  that  the  small  increase  in  budget 
was  sufficient  to  eliminate  almost  all  the  cases  of  infeasibility  (979  out  of  the  1000  instances  are  now 
feasible).  Thus,  we  conclude  that  the  violations  of  the  temporal  budget  constraints  were  of  small 
magnitude.  Also,  we  note  that  the  average  p  value  remains  rather  small  (0.0135).  This  exercise 
demonstrates  the  robustness  of  our  model. 


Value 

0 

.004 

.009 

.013 

.018 

.022 

.027 

.031 

.036 

Frequency 

303 

135 

114 

113 

68 

49 

47 

39 

22 

Value 

.040 

.045 

.049 

.054 

.058 

.062 

.067 

.071 

more 

Frequency 

16 

10 

11 

11 

8 

10 

8 

6 

9 

Table  4:  p  values  with  5%  increase  of  the  violated  temporal  budgets  constraints 
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Histogram  of 

(PertBasisOpt-  PertOpt)/  PertBasisOpt 


Figure  5:  Histogram  of  the  data  in  Table  2 


Count 

Average 

STD 

Max 

Min 

979 

0.0135 

0.0356 

1 

0 

Table  5:  p  values  with  5%  percent  increase  of  violated  temporal  budget  constraints 


Histogram  of 

(PertBasisOpt-  PertOpt)/  PertBasisOpt 


Bin 


Figure  6:  Histogram  of  the  data  in  Table  4 
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6  Concluding  Remarks  and  Extensions 


This  paper  addresses  a  resource  allocation  problem  faced  by  a  defense  agency  charged  with  devel¬ 
oping  CMs  in  a  dynamic  arms  race  against  an  adversary  that  seeks  to  cause  as  much  damage  as 
possible.  Such  scenarios  are  becoming  more  and  more  relevant  in  asymmetric  wars  between  gov¬ 
ernment  forces  and  insurgents.  The  main  contribution  in  this  paper  is  the  formulation  of  tractable 
network  optimization  models  that  encompass  the  essential  elements  of  the  problem.  The  models 
we  develop  are  deterministic  and  as  such  they  can  be  criticized  for  failing  to  address  the  (obvious) 
uncertainty  that  exist  in  the  development  times  of  the  CMs.  But,  the  extensive  numerical  analysis 
that  is  presented  in  Section  5  demonstrates  the  robustness  of  the  models  whose  outcomes  remain 
stable  when  there  is  some  “noise”  in  the  data. 

The  methods  we  developed  apply  to  a  modification  of  our  model  where  the  data  consists  of 
damage  rates  d^  for  rn  £  Jv[  and  W  C  W  and  the  damage  rate  in  the  presence  of  a  set  of  weapons 
W  and  a  set  of  CMs  M  is  =  minme»  d^  (in  the  presence  of  consistency,  (12)  shows  that  the 
model  we  study  is  an  instance  of  the  above).  When  all  weapons  are  operational  at  time  0,  this 
model  reduces  to  the  one  studied  in  Subsection  3.1  by  looking  at  W  as  a  single  weapon.  When  the 
weapons  are  not  necessarily  operational  at  time  0  and  weapon  w  becomes  available  at  time  sw  >  0, 
with  0  <  si  <  . . .  <  S|w|,  we  can  impose  a  modified  consistency  assumption  which  asserts  that 
the  ranking  of  the  CMs  against  the  sets  of  the  form  {1, . . .  ,w},  w  €  W,  is  the  same.  Under  this 
assumption,  the  analysis  of  Subsection  3.2  applies  and  the  decision  problem  can  be  reduced  to  (9). 

In  future  research  we  intend  to  explore  dynamic  versions  of  the  models  developed  here.  In 
particular,  we  intend  to  look  at  multi-stage  models  with  recourse.  That  is,  in  each  period  B  will 
be  able  to  observe  new  data  that  was  realized  since  his  previous  decisions  were  made  and  adjust 
the  decisions  accordingly.  We  plan  to  investigate  such  dynamic  models  in  both  a  “rolling”  and 
“folding”  horizon  frameworks  (in  the  rolling  horizon  framework,  each  decision  epoch  covers  a  fixed 
number  of  periods  in  the  future  while  in  the  folding  horizon  framework,  we  advance  towards  a  given 
target  date  and  so  the  decision  epochs  correspond  to  an  ever-decreasing  set  of  periods). 
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Appendix 


Proof  of  Lemma  2 

The  decision  version  of  the  constrained  shortest  path  problem  for  the  FCDS  graphs  defined  in  this 
paper  can  be  stated  as  following:  Given  positive  numbers  R  and  C,  determine  whether  there  is 
an  (O--D)-path  in  the  FCDS  graph  with  d-length  <  R  and  c-length  <  C.  This  decision  problem 
is  clearly  in  NP  since  the  length  of  a  path  in  a  graph  can  be  determined  in  polynomial  time.  We 
next  prove  that  this  decision  problem  is  NP-complete  by  showing  that  any  instance  of  the  Partition 
Problem,  known  to  be  NP-complete  (see  [6]),  can  be  reduced  in  polynomial  time  to  a  constrained 
shortest  path  problem  on  a  FCDS  graph.  Our  proof  modifies  arguments  of  [18]. 

The  data  for  an  instance  of  the  Partition  Problem  consists  of  a  set  of  n  >  1  positive  integers 
{ai,  02,  •  •  • ,  an}.  With  S  =  Y^i= l  aii  problem  is  then  to  determine  whether  there  exists  a  subset 

I  C  N  =  {1, . . . ,  n}  such  that  Yliei  a*  =  §■  We  next  construct  an  instance  of  the  CM  development 
problem  such  that  the  given  instance  of  the  partition  problem  has  a  solution  if  and  only  if  the 
corresponding  FCDS  graph  has  an  (O-D)-path  of  d-length  <  j  +  1  and  c-length  <  j. 

Without  loss  of  generality  assume  that  ai  >  02  >  03  >  •  •  •  >  an  >  1  (the  Oj’s  can  be  ordered 
using  0[n(lgn)]  comparisons).  Also  assume  that  S  is  even  and  ai  >  1  (else,  the  given  instance  of 
the  partition  problem  is  trivial),  implying  that  S  >  2.  Let  7  =  ^  <  1  <  S.  The  CM  development 
problem  that  we  construct  has  M.  =  {1, . . . ,  2n  —  1},  W  =  {1},  si  =  0  and  1C  =  {1}.  Hence,  the 
dependence  of  d{n  and  on  k  can  be  suppressed  and  these  are  given  by 

(Xl=\  Su  +  7,  g(gX7) ,  ^)  if  rn=2 i  -  1  for  i  G  {1, . . . ,  n  -  1} 

(35)  {tm,dm,cm)  =  <  (XX=i  Su,  0)  if  m=2i  for  i  G  {1, . . . ,  n- 1} 

(X-a=i  Su  +  7)  0,  an)  if  rn=2n  -  1 

(where  Su  stands  for  the  rt-power  of  S).  Also,  do  =  ^  and  T  =  X”=i  Su. 

Clearly,  0  <  ti  <  ...t2n- 1  <  T.  Also,  for  i  =  l,...,n-  1,  d2i- 1  =  sffib; j  >  |Itt  =  d2i 
and,  as  S{Sl  -  7)  -  Sl  =  Sl(S  -  1)  -  ai  >  0,  d2(i- 1)  =  §i  >  siS^)  =  Thus,  d0  >  di  > 

...  >  d2n- 2  >  0  =  d2n—\.  Following  (3),  the  corresponding  FCDS  graph  is  a  complete  directed 
graph,  i.e. ,  with  O  and  D  associated  with  0  and  2 n,  respectively,  its  vertex  set  is  V  =  {0, 1, . . .  ,  2n} 
and  its  edge  set  is  {( m,m ')  €  V  x  V  :  m  <  m'}.  Further,  for  each  edge  e  =  (m,  m')  for  which 
{is  G  V  :  m  <  is  <  rri'}  contains  an  even  integer,  say  £,  =  cm>  =  0  +  cm>  =  C(m  o  +  Cf£  m') 

and,  by  Lemma  1,  d(TO)7n')  >  d(m  o  +  dwmiy  Consequently,  such  edges  need  not  be  considered  in 
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exploring  the  existence  of  an  (O-D)-path  of  c- length  <  j  and  d- length  <  j  + 1.  Therefore,  one  can 
restrict  attention  to  paths  that  contain  all  of  the  even  vertices;  such  paths  are  determined  by  the 
set  of  odd  vertices  that  they  contain.  An  illustration  of  a  FCDS  graph  without  the  unnecessary 
edges  is  given  in  Figure  7. 


Figure  7:  The  FCDS  graph  constructed  for  the  proof  of  Lemma  2 

Consider  an  (O-D)-path  a  that  contains  all  even  vertices  and  whose  odd  vertices  consist  of 
{2 i  —  1  :  i  €  1}  where  /  C  J\f.  Then  c(a)  =  ai  and 


(36) 


«  =  e  <!>*+»£*+  e  yivy 

ieA/y  *  iei  b  iei\{n}  ^  7') 

=  (S  +  J2 


1(^-7) 


'Sl  Z— /  v  g 

iEl  iEl  i£l\{n} 


In  particular,  c{a)  <  if  and  only  if  <  I  and  d(a)  <  §  +  1  if  and  only  if 


(37) 


5 


|-1+E(|h+  E 

i£l  i£f\{n}  i£l 


If  I  =  0,  then  d(a)  =  S  >  §  +  1,  implying  that  this  case  can  be  excluded.  With  /  =7  0, 


(38)  0  <  E<^+  E  <ir><( 

i£l  i£l\{n } 


ai+l  \  ^  (  ^2i£AT  ai  \  ,Z)ieA/'\{l}  aC  _  <5  —  01 


s 


:)t  +  ( 


5 


)  =  7  + 


S 


=  1. 


Let  A  =  Ylieii §07  +  X)ie/\{n}(Or')’  As  the  ai’s  and  f  are  integers  and  0  <  A  <  1  (the  latter  by 
(38)),  f-l  + A  —  ai  ^  and  only  if  |r  <  °i>  he.,  (37)  is  equivalent  to  j  <  So, 
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c(cr)  <  ^  together  with  d(a)  <  §  +  1  are  equivalent  to  Yliei  ai  =  §•  Thus,  the  given  instance  of 
the  partition  problem  was  reduced  to  an  instance  of  a  decision  problem  of  a  constrained  shortest 
path  in  an  FCDS  graph.  □ 


The  general  model  used  for  the  numerical  experiments 


The  formulation  of  the  general  model  we  used  for  the  numerical  experiments  is  an  extension  of 
the  model  presented  in  Subsection  4.4  according  to  Subsection  3.3;  that  is,  the  time  at  which  the 
development  of  each  CM  starts  is  a  decision  variable.  Recall  that  in  this  model  we  refer  to  plans, 
where  each  plan  consists  of  a  pair  ( k,r ),  where  k  is  the  intensity  of  developing  the  CM  and  r  is 
the  start  time.  We  let  T  be  a  finite  set  of  all  possible  start  times.  We  assume  that  neither  the 
effectiveness  of  a  CM,  nor  the  cost  and  duration  of  developing  it  are  affected  by  r. 

As  in  Subsection  3.3  a  policy  ir  is  now  defined  as  a  set  of  triplets  (to,  k,  r)  with  distinct  values 
of  to  and  effective  policies  are  defined  in  terms  of  the  tm’^  (the  completion  times).  Consider 
the  modification  of  the  multiple-FCDS  graph  to  multiple-  T- FCDS  graph  where  vertices  and  edges 
are  defined  in  terms  of  the  triplets  (to,  k,  r)  instead  of  the  pairs  ( m,k ).  As  in  Subsection  3.3, 
if  the  development  plan  of  CMm  is  ( k,r )  then,  during  the  development  period  [r,  a  cost 

Cm’T^  =  C(k.T)  Per  unit  time  is  incurred  4.  The  total  expenditure  on  CMm  during  the  time  interval 

tm 

Ih  =  [T_hiTh]  is  c^’r)[min{ ~  max{Tft, r}]+.  For  each  edge  e  of  the  multiple-T-FCDS 
graph  and  h  =  1, . . . ,  H,  let 


(39) 


CmT)[mm{Th,t.£'T)  }-  rnax{T/j,  r}]+  if  e  terminates  at  (in,  k,  r) 
0  if  e  terminates  at  D. 


The  total  cost  associated  with  effective  policy  n  during  the  time-interval  Ih  is  then  expressed  by 
YleeE  cexe  aRd  this  sum  is  subject  to  the  corresponding  temporal  budget  constraint  Ch. 

Now,  consider  the  total  damage  associated  with  an  effective  policy  n.  Suppose  each  weapon  w 
becomes  operational  at  time  sw  >  0.  For  each  edge  e  of  the  multiple-T-FCDS  graph,  let  de  be  given 
by  (15  ),  where  the  d ™  values  are  given  by  the  variant  of  (14)  in  which  k  is  replaced  by  (k,r).  The 
total  damage  associated  with  policy  n  is  then  expressed  by  YIwgw  Yle£E(aw)  de-  Our  formulation 
is  given  as  follows: 

4In  our  experiments,  we  took  tm’T)  =  and  Cm’T ^  =  c^Vr,  i.e. ,  the  cost  rate  cS’r|  holds  during  the  development 

period  but  its  value  does  not  depend  on  the  start  time  r. 


28 


^2e<-(m,k,T)w  Xe  =  Y^e^(m,k,r)^  Xe  for  each  m  £  M,k  £  1C,T  £  T  and  W  £  W 

(40)  *e  =  1  =  Yje^D™  Xe  for  each  W 

xe  £  {0, 1}  for  each  e  €  E. 

In  particular,  if  configuration  a  =  {aw  :  w  £  W}  corresponds  to  x  satisfying  (40),  then 

(41)  Y  de  =  Y deXe' 

U)€W  e£E(crw)  e£E 

Given  such  a  solution  x ,  let  y(x)  =  [y{x)m'T\(m,k,T)&u  be  defined  by 

y(x)mT)  =  max  Y'  xe  for  each  (to,  A,  r)  £  U 

w  J 

e^(m,k,r)w 

(where  U  is  redefined  to  accommodate  triplets  (to,  k,  r)  6  M.  x  K,  x  T  instead  of  pairs  (to,  A)  € 
A4  X  JC);  equivalently,  y(x)  is  the  unique  solution  of 


y(x)mT)  >  Ee-(m,fc,r)« 

for  each  (to,  A;,  r) 

£  U  and  w  £  W,  r  €  T 

(42) 

y{x)m'  ■*  <  Y1W£V\>  J2e^(m,k,r)w  Xe 

for  each  (to,  A,  r) 

£  U 

y(x)mr)  e  {0,  1} 

for  each  (to,  A,  r) 

£  17. 

Now,  the  budget  constraints  are; 

(43)  Y  cexe  <  Ch  for  h  =  1, . . . ,  H 

e€E 

and 

(44)  XI  <  C. 

(m,k)eU 

The  problem  we  solved  in  our  numerical  experiments  is: 


(45) 


min 

x,  y 

s.t. 


Ydexe 

e£E 

(x,y)  satisfies  (40),  (42) (43),  and  (44). 
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